Galaxy SED catalogs

For select redshifts of \(z = \{6, 7, 8, 9, 10\}\) we produce galaxy spectral energy distribution (SED) catalogs for the flagship THESAN-1 simulation, directly mirroring the SUBFIND subhalo catalogs. The post-processing methodology for generating mock SEDs of THESAN galaxies from the star and dust distributions is outlined in Section 2.2 of Kannan et al. (2022).

Briefly, the intrinsic nebular emission is taken from Byler et al. (2017), which is based on photoionization calculations using CLOUDY. The underlying stellar radiation intensities and spectra are taken from the Flexible Stellar Population Synthesis (FSPS) model (Conroy, Gunn & White 2009). Dust attenuation is modeled using the Monte Carlo radiative transfer code SKIRT (Camps & Baes 2020), which accurately captures scattering and absorption by dust throughout the entire wavelength range considered 0.05 to 200 μm, discretized into 657 unequal bins. The radiative transfer calculations are performed on a Voronoi grid constructed from positions of the gas cells in the galaxy, mirroring the structure and geometry of the AREPO data as closely as possible. A redshift-dependent dust-to-metal ratio is used and dust emission is modeled assuming thermal equilibrium with the local radiation field. A Draine & Li (2007) dust mixture of amorphous silicate and graphitic grains, including varying amounts of polycyclic aromatic hydrocarbons (PAHs) particles, is assumed. Only sufficiently well-resolved, star-forming galaxies are considered, to ensure that the output SEDs from the SKIRT calculations are reasonably converged. We consider a galaxy well resolved if the stellar mass within twice the stellar half-mass radius is greater than 50 times the baryonic mass resolution of the simulation. We also require that at least one star particle be less than 5 Myr old, which ensures that there is recent on-going star formation within the halo. Given the resolution of the simulation, this roughly translates to haloes with a minimum star formation rate of just below \(0.1\,{\rm M}_{\odot} / {\rm yr}\).